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Effect of Improving the Lattice Gauge Action on QCD Topology * 

J. Grandy'' and G. Kilcup,*" 

^L-170, Lawrence Livermore National Laboratory, P. O. Box 808, Livermore, CA 94551 
''Physics Department, 174 West 18th Avenue, Ohio State University, Columbus, OH 43210 

We use lattice topology as a laboratory to compare the Wilson action (WA) with the Symanzik-Weisz (SW) 
action constructed from a combination of (1 x 1) and (1 x 2) Wilson loops, and the estimate of the renormalization 
trajectory (RT)[^ from a renormalization group transformation (RGT) which also includes higher representations 
of the (1 X 1) loop. Topological charges are computed using the geometric (Liischer's) and plaquette methods on 
the uncooled lattice, and also by using cooling to remove ultraviolet artifacts. We show that as the action improves 
by approaching the RT, the topological charges for individual configurations computed using these three methods 
become more highly correlated, suggesting that artificial lattice renormalizations to the topological susceptibility 
can be suppressed by improving the action. 



1. Improved Actions 

In lattice QCD, one of the fundamental chal- 
lenges is to minimize the errors caused by 
discretizing space-time. This is accomplished 
through a combination of advances in computer 
technology, and advances in the formulation of 
methods to solve the problem computationally, 
including the development of improved numeri- 
cal algorithms. We show an example of how an 
improved gauge field action can be used to sup- 
press artificial lattice contributions to physical 
measurements . 

We consider gluon actions that are constructed 
in a gauge invariant fashion, from a combination 
of Casimir invariants of closed Wilson loops. In 
principle, a lattice action of this type can consist 
of an arbitrary sum of Wilson loops, but a trunca- 
tion to a small set of localized loops is necessary 
due to computational expense. We study actions 
constructed from (1x1) and (1x2) loops: 

5 = i^ixi^ReTrl¥(i'^i) 
+ i^ix2 ^ReTrVK(ix2) 

+ ifixi ^Re[^(TrM/(ixi))2-iTrM^(i><i)] 

+ ^8^^^ E[?TrM/(i-)|^-i] 
where the actions and coefficients, in order of 



increasing improvement (approximation to the 
renormalization trajectory) are given by: 
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To compare these three actions, we generate 
four ensembles: 



Action 

Wilson 

SW 

SW 

RGT 



The two SW ensembles allow us to study the ef- 
fect of increasing f3, we have used estimates of cor- 
responding Wilson action (3 from the deconfining 
phase transition temperature calculation by Cella 
et al. Ql . Since we have used a modest number of 
configurations in each case, we focus on the qual- 
itative comparison between Wilson and improved 
actions. Further calculations with a larger num- 
ber of lattices would be needed for quantitative 
studies, for example, to determine the consistency 
and scaling of xt/mp. 
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2. Topology: Comparing Actions 

Lattice topology provides a test case for com- 
paring various gauge field actions. There are 
several prescriptions for measuring topological 
charge 

Q=^J d'xFix)F{x) 

on the lattice, and each prescription is subject to 
a different set of lattice cutoffs and renormaliza- 
tions which affect the measurement of the topo- 
logical susceptibility xt — (Q"^) tJ^6 pla- 
quette method the topological density F(x)F{x) 
is constructed from a product of lattice (1x1) 
Wilson loops. This method in general gives non- 
integer values of the topological charge, and is af- 
fected by large multiplicative and additive lattice 
renormalizations||]. The geometric method Q 
does guarantee an integer topological charge (ex- 
cept for "exceptional" configurations) but is not 
guaranteed to obey physical scaling in the contin- 
uum limit, and is in fact known to violate scaling 
for the Wilson action[^. Low-action dislocations 
which can be suppressed by improving the ac- 
tion]^ contaminate the geometric Xt- 

In the cooling prescription, ultraviolet fluctu- 
ations in the fields are removed by locally min- 
imizing the action in successive sweeps, isolat- 
ing instanton-likc configurations. After cooling, 
a single instanton configuration spanning several 
lattice spacings has a computed charge of nearly 
one using either the geometric or plaquette for- 
mula; we therefore apply the plaquette formula to 
the cooled configurations to obtain a value for Q. 
Lattice artifacts are very different among these 
methods, and we can in general get different re- 
sults for plaquette {Qp), the geometric (Qg), and 
the cooling (Qc) topological charges computed on 
the same original configuration. For improved 
actions, we expect lattice artifacts such as dis- 
locations to be suppressed, therefore we test this 
prediction by comparing the different topological 
charge methods with each other. 

The cooling prescription actually encompasses 
a family of cooling algorithms. Typically one 
cools by selecting a link U to minimize some ac- 
tion Sc , and since cooling is merely used as a tool 
to isolate instantons, there is no reason to tie Sc 



to the Monte Carlo gauge action S. The cooling 
algorithms Sc we consider here are linear com- 
binations of Wilson loops with coefficients C(ixi) 
and C(xx2)7 E^^d since action is minimized only 
the ratio ri2 = C(ix2)/c(ixi) is significant. The 
cooling algorithm with ri2 — —0.05 removes the 
leading scaling violation from the classical instan- 
ton action, and we also include cooling algorithms 
with ri2 = and ri2 = —0.093, which has been 
derived from a linear weak coupling approxima- 
tion to the RGT action, for comparison. For the 
case ri2 — 0, the lack of barrier to a decrease 
in the instanton size causes the instanton to dis- 
appear by implosion during the cooling process, 
and for ri2 = —0.093 a large instanton expands 
until halted by the boundary. We cool for 200 
sweeps for all three algorithms, and the compar- 
ison between these three in an indication of the 
systematic effect of picking some particular means 
of cooling. We note that with 200 sweeps of Wil- 
son cooling most of the topological charges are 
retained, since the large instantons haven't had 
enough time to implode. In general, we do not 
see any effect from the selection of the cooling 
algorithm, except perhaps in one ensemble. 

3. Results 

As described above, we compute Qp, Qg, and 
Qc on all of our lattices. We show two scatter 
plots (Figures |l], ||) highlighting the discrepancy 
between Qg and Qc- The best fit line is con- 
structed through the points on a scatter plot. The 
slope of this line is an estimate of the ratio of mul- 
tiplicative renormalizations, and should be close 
to 1 since both the geometric and cooling meth- 
ods give integer charges. The correlation 

_ {{Qg-Qg) (Qc-Qc)) 

Zgc — I 

\j{{Qg-Qaf){{Qc-Qc)^) 

between Qg and Qc is a measure of random ad- 
ditive artifacts seen by one method but not the 
other, such as dislocations which disappear dur- 
ing cooling therefore contributing only to Qg but 
not to Qc- The scatter plots show a strong cor- 
relation between Qg and Qc for the RGT action. 
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Figure 1. Comparison of geometric and cooled 
topological charges for Wilson action lattices at 
a ~ O.lfm. Each point represents one configura- 
tion, with the cooled charge as the absicca and the 
geometric charge on the uncooled lattice as the 
ordinate. The least squares linear fit is shown. 
Due to close overlaps there appear to be fewer 
that 35 points. 



and a far weaker correlation for the WA, suggest- 
ing that the effect of lattice artifacts on topologi- 
cal charge is far less pronounced for the improved 
action. 

In Fig. H we show a comparison between the 
WA, SW, and RGT actions of the correlation be- 
tween Qg and Qc- Fig. |^ similarly shows the 
correlation between Qp and Qc, computed in the 
same manner as Zgc, and numerical values are in 
Table |. The SW action serves as an intermedi- 
ate point between the other two actions, since the 
RGT action represents a better estimate of the 



Figure 2. Comparison of geometric and cooled 
topological charges for RGT action lattices at a w 
O.lfm. 



renormalization trajectory than the WA and SW 
actions. It is unclear whether the spread in Zgc 
at /3 = 4.43 is due to any systematic effect of the 
cooling algorithm. It is possible that for better 
improved actions, where the exponential suppres- 
sion of dislocations is greater than for the S-W 
action, increasing beta will have a more profound 
effect than we have seen here for the S-W action. 
For the plaquette method, we have shown previ- 
ouslyPl that the multiplicative Zp becomes less 



Table 1 

Correlations, SW Cooling Linear Fits 
Corr. WA SW, 4.2 SW, 4.43 RGT 
Zgc 0.28(13) 0.77(8) 0.80(6) 0.88(5) 
zpc 0.04(16) 0.29(16) 0.31(13) 0.32(16) 
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Figure 3. Correlation Zgc between geometric and 
cooled topological charges, for WA, SW, and 
RGT ensembles. Results for three different cool- 
ing methods are shown. 



severe as the action improves, and the increased 
correlation Zpc suggests that the additive renor- 
malization also decreases. In addition, improving 
the action is far more effective than increasing P 
for suppressing lattice artifacts. 

Future calculations can include a more com- 
prehensive study of other improved actions. 
Other methods for xt, including the fermionic 
method 0, and an indirect measurement by cal- 
culating the Tj' mass[H,||, should also be tested. 
Having established a correlation between cooled 
and uncooled topology and located individual 
instantons, we are now prepared to investigate 
Shuryak's picture of a dilute instanton gas, and 
the influence of instantons on hadronic physics by 
working directly on uncooled lattices. 
Acknowledgement. The calculations were per- 



Figure 4. Statistical correlation Zpc between pla- 
quette and cooled topological charges, for same 
ensembles, three cooling methods. 
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